Angle-resolved Landau spectrum of electrons and holes in bismuth 
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In elemental bismuth, emptying the low-index Landau levels is accompanied by giant Nernst 
quantum oscillations. The Nernst response sharply peaks each time a Landau level intersects the 
chemical potential. By studying the evolution of these peaks when the field rotates in three per- 
pendicular planes defined by three high-symmetry axes, we have mapped the angle-resolved Landau 
spectrum of the system up to 12 T. A theoretical model treating electrons at L point with an ex- 
tended Dirac Hamiltonian is confronted with the experimentally-resolved spectrum. We obtain a set 
of theoretical parameters yielding a good but imperfect agreement between theory and experiment 
for all orientations of the magnetic field in space. The results confirm the relevance of the Dirac 
spectrum to the electron pockets and settle the longstanding uncertainty about the magnitude of 
the g-factor for holes. According to our analysis, a magnetic field exceeding 2.5 T applied along 
the bisectrix axis puts all carriers of the three electron pockets in their lowest (j = 0) spin-polarized 
Landau level. On top of this complex angle-dependent spectrum, experiment detects additional and 
unexpected Nernst peaks of unidentified origin. 

PACS numbers: 71.70.Di, 71.70.Ej, 72.15.Gd 



I. INTRODUCTION 

Bismuth crystallizes in a rhombohedral structure. 
This is an intrigtiing alternative to the more symmet- 
ric cubic symmetry, favored by the balance in electron 
energeticsA^^-'l The departure from higher symmetry has 
made it a semi-metal with a significant role in the his- 
tory of metals physics (for reviews see^"^. In particular, 
it was in bismuth that a Fermi surface was experimen- 
tally determined for the first time seventy years agcP. 
This historical priority is not a consequence of intrinsic 
simplicity. The structure of Fermi surface in bismuth is 
quite unlike the sphere seen in the opening chapters of 
a standard condensed-matter textbook on the free elec- 
tron gas. On the other hand, and in spite of its complex- 
ity, this Fermi surface is one of the easiest to detect by 
experiment. Indeed, because of the lightness and high 
mobility of the carriers, quantum oscillations in bismuth 
are observable at relatively high temperatures and low 
magnetic fields. 

During several decades, numerous investigations of the 
de Haas- van Alphen and Shubnikov- de Haas effects were 
carried on bismuth in order to pin down the fine details of 
the Fermi surface2"i3. The well-established structure is 
sketched in Fig.l. It consists of one hole ellipsoid aligned 
perpendicular to the plane in which lay obliquely three 
slightly tilted-iS electron ellipsoids. Thanks to these in- 
tensive studies, the size of each ellipsoid and the tilt an- 
gle of the electron ellipsoids are known with a remarkable 
precisiorfiSl. 

In parallel to these experiments, a number of band cal- 
culations have been reported^^^^. However, the relevant 
energy gap of the system is in the range of 10 — 20meV 
and therefore beneath the typical energy resolution of 




FIG. 1: a) Sketch of the Brillouin zone and the Fermi surface 
of bismuth. The size of the Fermi surface is enhanced to make 
it visible. It consists of one hole pocket at the T point and 
three electron pockets at the L points of the Brillouin zone, 
b) In our experiment, 9 is the angle between the magnetic 
field and the trigonal axis, (f) is the angle between the binary 
axis and the projection of the magnetic field in the equatorial 
plane (defined by the binary and bisectrix axes) The three 
electron pockets are noted ei, e2 and 63 as seen in the figure 
according to their position to the angle (j>- 



ah initio calculations. The most successful band pic- 
ture has been a tight-binding model using semi-empirical 
parameterfPH. The results of this model are in rather 
good agreement with the details of the Fermi surface ac- 
cording to the experiment. 

Recently, there has been a renewal of interest in the 
electronic spectrum of bismuth, triggered by t he rela- 
tively easy access it provides to the quant um limitP^ISSIll] 
and the presence of Dirac fermions^Si^fil. The quantum 
limit is attained when the magnetic field is strong enough 
to confine electrons to their lowest Landau level. For 
ordinary bulk metals, the size of the required magnetic 
field is well beyond the limits of current technology. In 
a low-density semi-metal such as bismuth, on the other 
hand, a modest magnetic field applied along the high- 



symmetry axis known as the trigonal axis allows one to 
attain this limit. High- field exp eriments reporting peaks 
in the Nernst coefficieniFSElIll]^ jumps in magnetization^^ 
and Hall plateauP^ raised the issue of collective effects 
in a three-dimensional electron gas at high magnetic 
fieldpHal. 

In the context of these unexpected experimental find- 
ings, the high-field phase diagram o f bism uth became 
subject to a theoretical reinvestigatiorp2E31. The Landau 
level spectrum of a Fermi liquid is determined by the de- 
tailed topology of the Fermi surface. In bismuth, this 
spectrum becomes complex in the vicinity of the quan- 
tum limit for several reasons. To preserve charge neutral- 
ity, any modification in density of one type of carriers, 
caused by emptying a Landau level, should be compen- 
sated with an equal change in the density of the carriers 
of the other sign. The chemical potential shifts as the 
magnetic field is scanned and quantum oscillations are no 
longer a periodical function of B^^. Moreover, the large 
spin-orbit coupling leads to an angle- dependent Zeeman 
energy. For both electrons and hole, the anisotropy of 
the "spin mass" is different from the anisotropy of the 
cyclotron mass. As early as 1964, Smith et al. investi- 
gated this feature^. 

Dirac fermions provide a major motivation to revisit 
the Landau spectrum of bismuth. It has been known 
for a long-time that the energy-momentum dispersion of 
electrons at the L-point is not parabolic and that Dirac 
Hamiltonian is appropriate to treat therrpSHSS] Dirac 
fermions have been identified as the source of the en- 
hanced diamagnetism in bismuttpSI^ and may generate an 
unusual weak-field Hall conductivity^. During the past 
decades, several theoretical models of three-dimensional 
Dirac fermions have been p roposed and confronted with 
the experimental datgP^^I^SHlI]^ i^^t in contrast to the well- 
established topology of the Fermi surface^°, the details 
of the relevant Hamiltonian has not been pinned down. 

During the last few years, bismuth has attracted at- 
tention from other angles of research. In particular, in- 
frared conductivity measurements have detected a collec- 
tive electron-plasmon mode (a "plasmaron" pSl and the 
insulating-like magnetoresistance has been a subject of 
debate and invest igatioiP^ES. 

In this paper, we report on a study of angle-dependent 
Nernst effect in bismuth, which maps the evolution of 
Landau levels of electrons and holes in the whole solid 
angle up to 12 T. The data are compared to the results 
of a theoretical calculation based on a model with an ex- 
tended Dirac Hamiltonian. The comparison between the- 
ory and experiment narrows down the choice of parame- 
ters for this model. Two previous reports have mapped 
the angle-dependent Landau spectrum of bismuth in a 
restricted angular window. Smith and co-workers^ used 
Shubnikov- de Haas measurements to determine the Lan- 
dau spectrum in (trigonal, binary) and (trigonal, bisec- 
trix) planes. This pioneer work established the bulk of 
electronic properties, but left unsettled the precise mag- 
nitude of the hole Zeeman energy. Kajimura and co- 



workers^ used giant oscillations of ultrasonic attenua- 
tion to obtain a high-resolution map of the Landau spec- 
trum near the trigonal axis. Our study, continuing this 
line of exploration, uses giant quantum oscillations of the 
Nernst coefhcient to obtain a high-resolution map in the 
whole solid angle. The verification of the non-interacting 
spectrum in bismuth provides a basis to find an explana- 
tion for the additional features (i.e. Nernst peaks^^ ^^ ^51) 
which are not expected in this picture. 

Among the details of the one-particle picture coming 
out of our investigation, let us highlight two. First, our 
results definitely settle the uncertainty on the magnitude 
and anisotropy of the Zeeman energy for holes. Second, 
when the field is along the bisectrix, the spectrum of the 
three electron pockets is very close to a perfect Dirac 
spectrum where the cyclotron and Zeeman energies are 
indistinguishable. In other words, spin and orbital angu- 
lar momenta are described by a single quantum number. 
In this configuration, a magnetic field as low as 2.5 T is 
enough to attain the quantum limit of all electron pockets 
to put them in their lowest spin-polarized Landau sub- 
level. In other words, all electron-like carriers become 
spin polarized, residing in their lowest {j = 0; where j is 
the total angular momentum) Landau level. 

In addition to the case of bulk semi-metallic bis- 
muth, the results reported here may have implications 
for the analysis of the electronic properties of topologi- 
cal insulators^. In particular, in the case of Bii_a;Sba; 
alloys , the electronic spectrum of bulk bismuth and 
bulk antimony is the starting point to analyze the surface 
states^ of these systems. 



II. THEORETICAL TREATMENT OF THE 
SPECTRUM AT L POINT 

The essential properties of electrons at L point are 
the following: (1) The Fermi surface is ellipsoidal and 
tilted. (2) The electronic dispersion is non-parabolic. 
(3) The Landau levels display a slight spin splitting. 
The first property was already realized at the early stage 
of investigatiorP. The second one, the non-parabolicity 
of the dispersion, can be naturally obtained on the basis 
of the k ■ p perturbation theory. The k ■ p method 
was first applied to bismuth by Cohen and BlounPS. 
Soon afterwards, Lax and co-workers used this model 
to analyze their experimental datcpS. Because of 
the narrow gap character of L point, the number of 
bands needed is small and a model with coupled two 
bands gives a very good agreement with experiment. 
The k ■ p theory for a two-band model without spin 
gives the energy under a magnetic field in the form: 

E = ±[A^ + 2A{{n+l/2)nuj^ + h^kl/2m,}Y^^ , 
where ujc is the cyclotron frequency and rriz is the 
effective mass along the direction o f the m agnetic field. 
The third property, the spin splitting^SSEI]^ pan be taken 
into account if we introduce the spins with an effective 
g-factor (Fig. [2] (b)). The energy now becomes E = 



(a) k.p 
w/o spin 



(b) SBR (c) Dirac 



(d) Extended 
Dirac 




FIG. 2: Evolution of the bismuth model. The horizontal lines 
indicates the position of Landau levels with opposite spins. 



± [A^ + 2A{{n+ 1/2) hujc + h'^kl/2m^ ± g^fffiBB/2}] ^^^ , 
where /iB = \e\h/2mcC is the Bohr magneton. In spite 
of its very simple form, this two-band model, which we 
will call the Smith-Baraff-Rowell (SBR) model, gives 
a very good account of the experiment, in particular 
wh en the magnetic field is oriented close to the trigonal 
g^-^^ j8 | 34 | 3 6] However, there is a serious intrinsic problem 
in the SBR model. When the term —gcsiJ,BB/2 is 
large and negative, the energy becomes imaginary. As 
a consequence, it is difficult or impossible to use this 
model when the field is large and oriented far off the 
trigonal axis. 

In the original Cohen-Blount theory, the k ■ p theory 
was derived including the spin-orbit coupling (not fully 
contained in the SBR model) , which is extremely large in 
bismuth {Eso - 1.5eV20). Wolff^^^ found that the Cohen- 
Blount theory is essentially identical to the Dirac theory 
in the relativistic quantum mechanics. This finding led 
to the birth of the Dirac electron in solids. The two band 
kp model with spin-orbit coupling (Dirac modeP^^) gives 
the energy^in terms of a quantum number j = n+l/2+s 

as: E = ± [A^ -|- 2A {jhuj^ + fi^kl/2m^]] ^''^ . With this 
Dirac model, the spin-orbit coupling is taken into account 
and the problem of imaginary energy is removed, so that 
it can be safely used for all field directions. On the other 
hand, in this model, the spin splitting is exactly the same 
as the Landau level splitting, and a double degeneracy 
occurs for all Landau levels except for the lowest Landau 
level j — 0(Fig. [21(c)). This is not in agreement with the 
experimental observation of spin splitting. 

To repair this discrepancy and to understand the spin 
splitting, Baraff extended the Dirac model by consider- 
ing the effect of other bands outside the closely coupled 
two bands based on the perturbation theorySH. As a re- 
sult, he showed that it is possible to produce a difference 
between the orbital level splitting and the spin splitting 
(Fig. [2] (d)). Although the Baraff 's model is complex to 
use for the interpretati on of exp erimental results, Dres- 
selhaus and co-workergSS^ESMl have succeeded to sim- 
plify the Baraff's model and to obtain theoretical results 
in good agreement with experiments at least when the 
magnetic field is oriented along the binary and bisec- 



trix axes. In the simplified model, the energy is given 

hy E = ±[h? + 2A [jhuj^ + h^Pj2m,}] ^'"^ ± g'/iBB/2, 
where the additional g-factor expresses the effect of out- 
side bands. With this extended Dirac model, all three 
properties above are taken into account, and no imagi- 
nary energy appears in the entire solid angle. The num- 
ber of parameters is the same as the SBR model and it 
is as easy to handle as the SBR model. Consequently, 
the extended Dirac model will be the most appropriate 
model for the present purpose. 

It should be noted here that the energy band model 
of bismuth has been the subject of an alternative ap- 
proach, which starts with a simple cubic lattice without 
spin-orbit coupling. The Hamiltonian for bismuth is then 
expanded to first order with (i) the distortion of the lat- 
tice and then with (ii) the spin-orbit coupling, and (iii) 
the distance in k space from the L point. This approach 
was first developed by Abrikosov and FalkovskjEl, and es- 
tablished by McClure^^. One important consequence of 
this approach is that the dispersion contains in addition 
to a term linear in fc term a quadratic term. This is the 
so-called non-ellipsoidal non-par aboli c dispersion equiva- 
lent to the form found by CoherP^ESI, xt is also the char- 
acteristic of this model that the Hamiltonian is directly 
related to the relative displacements of the atoms, the 
rhombohedral deformation, and the spin-orbit couplings, 
so that we can predict the effect of substitution of other 
elements or the effect of pressure. The obtained disper- 
sion is similar to the Dirac model except for the presence 
of quadratic (non-ellipsoidal) term. However, there is no 
compelling evidence to prefer the non-ellipsoidal model 
over the ellipsoidal modeP, in spite of various attempts 
to use the non-ellipsoidal model for the analysis of exper- 



imental dati 
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Moreover, it is not clearly understood 



whether this non-ellipsoidal model can reproduce the ex- 
perimentally resolved spin splitting. 



III. THE MODEL 

The effective Hamiltonian for electrons in bismuth at 
L point in the Brillouin zone can be described in terms of 
the two-band model (written as a 4 x 4 matrix including 
the spin degr ees of freedom). The Dirac Hamiltonian is 
given bjlaSEil 






A i 

KKct 



AK a- 

A 



(1) 



where cr is the Pauli matrix and the 2x2 unit matrix 
is omitted. When the direction of the magnetic field is 
chosen as direction z, and the directions x and y are 
chosen to be perpendicular to each other, the (quasi) 



TABLE I: All theoretical parameters for the electron mass 
rrii and the additional g- factor gl, and for the orbital- and 
spin mass of holes, which give the best fitting results with 
experiments. The tilt angle of the electron pocket is 6.2°. 



Electrons 


mi 


m2 


ms 


7714 




0.00124 


0.257 


0.00585 


-0.0277 




9i 


9'2 


9'3 


5i 




-7.26 


24.0 


-7.92 


9.20 




Holes 




Ml = M2 




Ms 


orbital mass 
spin mass 




0.0678 
0.0319 




0.721 
200 



gap at L point (2A) 
electron-hole hybridization (-Eo) 



15.3 meV 
38.5 meV 



TABLE IL Various physical values along the principal axes. 
nic and Mc refer the cyclotron mass of electrons and holes. 
niz and M^ are the band mass of electrons and holes along the 
field orientation. They are all normalized by the bare mass of 
electron, me. 







B 1 Binary 


B 1 Bisectrix 


B II Trigonal 


mt^ 




0.0272 


0.00189 


0.0125 


mf'^^ 




0.00218 


0.00375 


0.0125 


mf 




0.00124 


0.257 


0.00585 


mf'"'' 




0.193 


0.0653 


0.00585 


f' 




73.5 


1060 


159 


^el,e3 




917 


533 


159 


5'" 




-7.26 


24.0 


-7.92 


g/cl,o3 




16.2 


0.545 


-7.92 


1 + {m,g'Y 


V2 


0.90 


1.02 


0.950 


1 + {m^g'Y 


l,ed/2 


1.01 


1.00 


0.950 


Mc 




0.221 


0.221 


0.0678 


M, 




0.0678 


0.0678 


0.721 


G 




0.791 


0.791 


62.6 


Ez/Tiujc 




0.0875 


0.0875 


2.12 



wave vector K is related to the crystal wave vector by 



K^ ± iK„ 



K,= 



nkz 



(2) 
(3) 



with the cyclotron mass nic, the longitudinal mass m^, 
kx = kx + {eB/2c)y and ky = ky — {eB/2c)x, i.e., the 
symmetric gauge. The wave vector K satisfies the com- 
mutation rule: 



K X K = 



iheB 



rricC 



(4) 



In Wolff's version of Dirac theory, the eigenvalues of this 
Hamiltonian can be obtained by considering its squared 
equation. 



_2, /A2+A(K.<t)2 

^o^^[ A2 

The quantity {K ■ cr)'^ is written as 
he 



AiK-a)'^"^ 



{K-a) 



K' 



(T B, 



rricC 



(5) 



(6) 



by using the commutation rule. Then, the squared wave 
equation becomes 



f,^* 
I Jf ' 



v- 



E^ 



2A 



-V', 



where 



,5^* 



K^ 
^ 



gfiBS ■ B. 



Here we introduced the effective g-factor 
„rne 



9 



nir. 



(7) 



(8) 



(9) 



with the spin operator s = cr/2. If wc introduce a "non- 
relativistic" energy E' — E — A (^ A), the right hand 
side of Eq. ([t]) is approximated as 



A^ 



2A 



tP ~ E'^/j. 



(10) 



Therefore, the Hamiltonian Jf* is the effective Hamil- 
tonian of J^o in the non-relativistic limit. It should be 
emphasized that the g-factor of J^* is written by the cy- 
clotron mass rric, namely, the spin mass of J^* is exactly 
the same as its orbital mass. This is one of the most im- 
portant signatures of a Dirac spectrum. J^* is just the 
Hamiltonian of free electrons, so that its eigenvalue is : 



e{n,s,kz) = [n+ - \ htOr. 



h'kl 

2m-, 



sg^BB^ 



2 1,2 



SJI3B, 



h'k_ 

2m, 



(11) 



where /3 = huJc/Bz- Finally, the energy of J^q is given by 



E = ±y/A^ + 2Ae{n,s,k^), 



(12) 



where the "-f " and "— " signs correspond to the conduc- 
tion and valence band, respectively and the origin is set 
at the center of the gap at L point. 

The Hamiltonian J^ is obtained by restricting the 
number of the bands under consideration to two. How- 
ever, in addition to these two bands, there are other 
bands and the existence of these outside bands splits the 



degenerate Landau levels of the Dirac model. By means 
of the perturbation theory, this effec t can be described 
up to the first order of B bjSHIlIllIl 



-^P = g'fJ-BS ■ B. 



(13) 




Then, the total Hamiltonian of the extended Dirac model 
is given by J^f — J% + J^,, and the total energy is by 



En,: 



where /3o = \e\h/mcC = 0.1158 meV/T. Note that, for 
the lowest Landau level (j = 0), the effect of the out- 
side bands appears in more complex form as is discussed 
by Vecchi and co-workers analyzing optical conductiv- 
ity datcP2. However, in contrast to the optical measure- 
ments, DC transport does not probe the evolution of the 
lowest Landau level. Therefore, here we express the ef- 
fect of outside bands effect just by J^ for j > uni- 
formly. We also note that g' can be positive and negative 
as in the theory of BarafP^, and Maltz and co- workers^, 
while i t app ears as \g'\ in the theory of Vecchi and co- 
workerPSSSl^ ^j^g present experimental results show that 
the level cross between 0+ and 1^ cannot be interpreted 
by \g'\. The obtained Hamiltonian of the extended Dirac 
model has a quite similar form as the Hamiltonian used 
by Alicea and Balentd^^l, though the fitting parameters 
are different. 

When jmc/3QB/mcA <C 1, we can employ "non- 
relativistic approximation" as 



2 1,2 



En,: 



± 



2m, 



rric 
nir 



1 

+ 2 
PoB. 



-f^oB 



(15) 



Then the ratio of the Zeeman splitting to the orbital split- 
ting is estimated as 



Ez^ 

huic 



Er, 



£"»).- 



(TOe/mc + g'/2) I3qB 



E., 



■"+i,± 



E, 



n,± 



1 



{mc/mc)f3oB 

2?71„ 



(16) 



The values of 1 -I- rn^g' /2m. 
should be noted that Eq, 



listed in Table |l!) It 

( |16[ ) is obtained on the basis 

This is a good ap- 

Therefore, 



of the non-relativistic approximation 

proximation only when jmc/3oB/m,cA. <^ 1 

the estimation based on this approximation is only valid 

for relatively large rric, e. g. when B// trigonal. For 

other configurations its domain of validity is restricted 

to the low magnetic fields. 

If we consider the angular dependence of B, it is very 
useful to describe m.^ in terms of the mass tensor m as 



det m 



(17) 



in which the longitudinal effective mass m,z is found by 
taking the component of m along the magnetic field. 



b ■ m ■ b. 



(18) 



where b is a unit vector along the magnetic field. The 
effective mass tensor for electrons at L point is given by 



'mi 
m — rrin I m2 m^ 

7714 7713 , 



(19) 



where 1,2, and 3 refer to the binary, bisectrix, and trig- 
onal axes, respectively. The tilted angle 9 is written as 
tan 261 = 2to4/(to2 -I- TO3). The other two electron pock- 
ets are obtained by rotations of ±120° about the trigonal 
axis. The additional g-factor is given in the same manner 
as 



9 



bg'b, 

'g'l 0' 

9i 93J 



(20) 
(21) 







Now we discuss the energy dispersion of holes at T 
point. In a Face-Centered-Cubic (FCC) lattice struc- 
ture, the two T points are strictly equivalent to the six L 
points. It is natural, therefore, to consider the holes at 
T point within the same framework as the electrons at 
L point. Note, however, th at the b and gap at T point is 
very large Ay (^ 200 meVE^IHini) ^ which is an order of 
magnitude larger than the gap at L point. Therefore, the 
non-relativistic treatment of holes becomes a reasonable 
approximation. Then the hole energy can be expressed 
like J^* for the electrons. Adding the additional Zeeman 
term, which originates from the outside bands, we obtain 
for holes an expression similar to the one used in the SBR 
model^^'': 



En 



A- E = 



hQc 



2M, 2^^ ' 



(22) 



where Eq is the hybridization energy. The cyclotron fre- 
quency (ric), the longitudinal- (M^) and cyclotron-mass 
(Mc), and the g-factor (G) for holes is given in the same 
manner as for electrons. The Zeeman splitting energy of 
this model is: 

Ez = GfiBB, (23) 

and so the ratio to the orbital splitting is given by: 

Ez GM, 



huJr 



2mr 



(24) 



This time, this formula gives good approximation for 
whole direction and amplitude of B, since the holes are 



non-relativistic, unlike in the case of electrons, Eq. ( 16 1 



The effective mass tensor has the form 



M 



Ml 











A/2 











Ma 



(25) 



The effective g-factor is defined in terms of the spin mass 

Ms by 



G = 2m, 



b-Ms-b 

detMs 



(26) 



This effective g-factor already includes the effect of the 
outside bands, which is closer to the initial two hole bands 
than that for clcctrontP^. So, the effective g-factor of 
holes may be much larger than that of electrons. 

The Fermi level E-p is determined by the charge neu- 
trality condition: 



J2NedEF)^Nh{EF) 



(27) 



i=l 



where A^e,i and Nh are the carrier numbers of electrons 
at i-th pocket and holes, respectively. The mass and g- 
factor tensors used in the present work are listed in Ta- 
ble [l] The components of the mass tensor are comp arable 
but not identical to those previously reportecpISEl]. They 
have been tuned to give the best agreement with the ex- 
perimental data obtained in this study. The parameters 
A = 7.65 meV and Eg = 38.5 meV are the same used in 
the SBR modeP^. 



induced by magnetic field is negligible and the field de- 
pendence of the Nernst coefficient is set by the variation 
in the transverse electric field as the field is scanned in 
presence of a constant heat current. The transverse volt- 
age was measured by an EM electronics Nanovoltmeter 
and collected through copper wires with no junction on 
their way between room temperature and 4.2K. In this 
way, thermoelectric noise was suppressed to measure a 
DC voltage of the order of 1 nV. Measurements were 
done in a dilution refrigerator inserted in a 12 T super- 
conducting magnet. 

The relative orientation of the crystal and the magnetic 
field was controlled using two coupled Attocube rotators. 
A first rotator allowed to explore a window of ±60 degrees 
in the principal rotational plane. A second smaller rota- 
tor allowed to rotate in a plane perpendicular to the prin- 
cipal plane of rotation and was used to adjust the plane 
of rotation. Angles were determined by a pair of perpen- 
dicular Hall probes. In this study we exclusively focus on 
the angular evolution of the field position of the Nernst 
peaks in order to extract the angle-dependent spectrum. 
We have repeated the experiment on several crystals to 
check the reproducibility of the results. A Nernst volt- 
age is the vector product of the thermal gradient and the 
magnetic field. In these experiments, the heat current 
was applied along the rotation axis. Therefore, during 
the rotation, the thermal gradient was perpendicular to 
the magnetic field, but not the electric field measured 
by a set of two fixed electrodes. In order to scan more 
than 90 degrees, we measured the electric field along two 
perpendicular directions using two pairs of electrodes. 



IV. THE EXPERIMENTAL PROBE AND 
SET-UP 



In the vicinity of the quantum limit, Nernst effect be- 
comes extremely sensitivity to quantum oscillationa^. 
Giant oscillations of the Nernst response have been ob- 
served in three semi-metals: bismuth^^ ^'', Bio.geSbo.oJ^ 
and graphite^. Each time a Landau level intersects with 
the chemical potenti al, Ne rnst signal sharply peaks. Re- 
cent theoretical worka^SSIl have explained the large mag- 
nitude of the quantum oscillations in this context. 

The giant oscillations dominate a monotonous back- 
ground which is unusually large. Indeed, the magnitude 
of the low-field Nernst coefficient in bismuth peaks to the 
large value of 7meVK^^ at about 3 K^. This is three 
orders of magnitude larger than the typical vortex Nernst 
signal of a type II superconductor. Thus, a Nernst exper- 
iment on a bismuth single crystal consists in measuring 
a relatively large transverse voltage of a few microvolts 
generated by a temperature difference as small as few 
miliKelvins. 

We employed a standard one-heater-two thermometers 
setup to measure the Nernst effect. In our range of inves- 
tigation, phonon thermal conductivity exceeds the elec- 
tronic heat conductivity by several orders of magnitude. 
As a consequence, the change in thermal conductivity 



ALONG THREE HIGH-SYMMETRY AXES 



Theoretical results 



The bottom positions of Landau levels as a function 
of B^^ are plotted in Fig. [3J As the field is swept, the 
levels of electrons move upward and those of holes move 
downward. The Fermi level shifts with this change in 
the density of electrons and holes in order to satisfy the 
charge neutrality. 

As seen in different panels of Fig. [SJ for _B j] bisectrix 
and B \\ binary, the Fermi level moves downward with 
increasing i?, whereas it moves upward for B \\ trigonal 
as a consequence of the rising Oj^ level. 

The variation of carrier density with magnetic field is 
shown in Fig. HI The density of carriers basically in- 
creases with increasing magnetic field due to the Landau 
degeneracy. Only for B \\ binary, the carrier density of 
the electron pocket along the bisectrix axis decreases with 
increasing magnetic field since the level of 0~2 approaches 
to E-p as seen in Fig. Islc). 
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FIG. 4: Carrier density as a function of B ^ for a magnetic 
field along a) Trigonal, b) Bisectrix and c) Binary direction. 



FIG. 3: Landau levels (fc^ = 0) as a function of B~^ for a 
magnetic field along a) Trigonal, b) Bisectrix and c) Binary 
direction. The origin of the energy is taken at the bottom of 
the conduction band at zero field. 



B. Experimental results 

In Fig. [5j the quantum oscillations of the Nernst re- 
sponse are plotted against the inverse of magnetic field 
for three different axes. In the same figure the expected 
theoretical positions for the intersection of the Fermi level 
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FIG. 5: Comparison between theory (vertical lines) and ex- 
periment (solid curves) for a magnetic field along a) Trigonal, 
b) Bisectrix and c) Binary direction. For each configuration, 
the Nernst voltage is plotted as a function of B^^). 



with each Landau level is marked by vertical lines. 

When the field is along the trigonal axis (the top 
panel), experiment detects oscillations with a relatively 
simple pattern as previously reportecpS. The Nernst re- 
sponse is dominated by holes, which give rise to peaks 
with a periodicity of 0.147 T^^. At high enough field 



{B > 1.5T) (that is for low Landau-level numbers), elec- 
tron peaks become detectable and the Zeeman splitting 
of holes becomes apparent. The quantum limit of holes is 
marked by a sharp peak at 9 T (B~-'^=0.11 T). Above this 
field holes reside in their lowest spin-polarized Landau 
level and electrons are at their two lowest Zeeman-split 
Landau level. 

When the field is perpendicular to the trigonal axis, the 
Nernst response is no more dominated by holes. As seen 
in the middle panel Fig. [5) large Nernst peaks due to elec- 
tron pockets become visible when the field is along the 
bisectrix. In this configuration, a field as small as 2.5 T 
puts the carriers of all three electron pockets in their low- 
est spin-polarized Landau level (that is 0^). Only above 
this field, the hole peaks become prominent. When the 
field is along the binary axis (the bottom panel of Fig. 
[5]), the quantum limit of two out of the three electron 
pockets is attained at 1.5 T. Hole quantum oscillations 
dominate the spectrum, as the field exceeds this magni- 
tude. The quantum limit of the third electron pocket is 
attained in a field of about 11 T. 

A careful examination of Fig. [5] indicates that for both 
holes and electrons, the agreement between experiment 
and the theory is rather good. 

The B II bisectrix configuration deserves special atten- 
tion. The upper panel of Fig. [6] displays the Nernst 
voltage as a function of B^^ for this configuration in a 
limited field range (0.3T < B < 3T) in order to highlight 
the Dirac spectrum of the electrons. As seen in the figure, 
quantum oscillations of the electron pockets are strictly 
periodic as a consequence of the stability of the chemical 
potential in this field window. Nernst peaks appear with 
remarkable regularity at multiples of 0.405 ± 0.005 T~^. 
The simplicity of this spectrum is a result of two inde- 
pendent factors. First, in the B || bisectrix configuration, 
the cross-section of the two pockets, el and e3 is exactly 
twice (that is l/cos(7r/3)) the cross section of the pocket 
e2. Second, in this configuration, the Zeeman energy and 
the cyclotron energy are almost identical, in other words, 
the n,-f and n-j-l,- states become almost degenerate. 

As a consequence of these two features, Nernst peaks 
for e2 emerge at multiples of 0.81 T~^ and those for the 
pockets el and e3 at multiples of 0.405 T~^. Therefore, 
the experimentally resolved peaks at odd or even mul- 
tiples of 0.405T~-'^ are respectively fourfold and sixfold 
degenerate. 

As seen in the Table II, when B || bisectrix, the theoret- 
ical magnitude of 1-1-^'^'^ remains close to unity leading 
to a spectrum very close to a purely Dirac spectrum in 
agreement with the experiment. As seen in the lower 
panel of Fig. [g] when the B~^ position of the Landau 
levels are plotted vs. the index number, the intercept is 
very close to zero as expected in the case of a the Dirac 
spectrum. 

We note that in this configuration, when the field ex- 
ceeds 2.5 T, electron- like carriers will reside in their low- 
est spin-polarized Landau level. 
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FIG. 6; Top: Nernst voltage as a function of B^^ with the 
magnetic field along the bisectrix. Electron peaks occur with 
a periodicity of 0.405 T~^ . This period is twice the cross sec- 
tion of an electron ellipsoid along the bisectrix axis. When, 
B~^ is an even (odd) number of times this value, a Nernst 
peak is sixfold (fourfold) degenerate. Bottom, the B~^ po- 
sition of Nernst peaks as a function of their Landau index 
for the three pockets. Dirac spectrum leads to a vanishing 
intercept. 



VI. ANGLE-RESOLVED LANDAU SPECTRUM 

The angle-dependent Nernst effect was studied for a 
magnetic field rotating in three perpendicular planes of 
rotation. Fig. [7] presents a typical set of Nernst data. 
The upper panel (Fig. JTJa) presents the Nernst signal, 
Sa;j, as a function of magnetic field for different field orien- 
tations tilted off the trigonal axis from -40 degrees to 40 
degrees. Quantum oscillations are clearly visible. Each 
peak corresponds to the intersection of a Landau level of 
one of the pockets with the chemical potential. The posi- 
tion of peaks shifts as the magnetic field rotates. A color 
plot of the same data is presented in the lower panel and 
reveals the Landau spectrum. The angular evolution of 
the Nernst peaks plots bright lines which take the shape 
of a network. The quasi-horizontal lines which curve up- 
ward as the field is tilted off the trigonal axis correspond 
to the Landau levels of the hole pocket: as the field is 
tilted, the cross section of this ellipsoid increases. The 
quasi-perpendicular lines correspond to the three elec- 
tron pockets, which lay almost perpendicular to the hole 
pocket. 

We used the same procedure to determine the color 
map of the Nernst signal in three principal planes. The 
results are presented in Fig. [8| In each panel, there are 
two adjacent parts representing the data obtained from 
different pair of electrodes on the same sample. No dis- 



FIG. 7: a) The evolution of the field dependence of the Nernst 
signal Sxy with the orientation of the magnetic field as the 
field rotates in the (trigonal, binary) plane. Curves are shifted 
for clarity, b) Color map of the data presented in the upper 
panel. Bright lines track the angular evolution of the Nernst 
peaks corresponding to the intersection of a Landau level and 
the chemical potential. 




FIG. 8: Color map of S^,,, when the magnetic field rotates in 
three perpendicular planes: a) (Trigonal, Binary) plane (i.e. 
(j) = 0°), b) (Trigonal, Bisectrix) plane (i.e. (j) = 90°) and c) 
(Binary, Bisectrix) plane (i.e. 6 = 90°). On top of the color 
map, theoretical lines corresponding to the hole (purple) and 
the three electron pockets, ei (dotted), e2 (black solid) and 
63 (gray solid) are plotted. In each panel there are two parts 
corresponding to the data obtained with two different pair of 
electrodes. 
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FIG. 9: The Landau spectrum for a magnetic field rotating in 
the (trigonal, bisectrix) plane according to theory (lines) and 
experiment (symbols). Additional peaks refer to those which 
can be clearly identified as not associated with an electron 
Landau sub-level. 
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FIG. 10: The Landau spectrum for a magnetic field rotating 
in the (trigonal, binary) plane according to theory (lines) and 
experiment (symbols). Additional peaks refer to those which 
can be clearly identified as not associated with an electron 
Landau sub-level. 



continuity appears as one crosses from one set of elec- 
trodes to the other. Theoretical lines derived from the 
model discussed above are put on the top of the color 
plots. As seen in the figure, these lines are in rough 
agreement with the experimental data. 

The three following figures allow a comparison between 
the theory and experiment for each of the three rotating 
planes. In each figure, the angular variation of the field 
at which a Landau band crosses the chemical potential is 
compared to the angular variation of the fields at which 
the Nernst signal peaks. 

It is easy to distinguish the Nernst peaks associated 
with holes, since there is a single hole pocket and their 
intensity evolves continuously with rotation. As seen in 
the figure, the theoretical model gives a very good ac- 
count of the experimentally-resolved hole peaks in all of 
the three planes. When the field rotates in the (binary, 
bisectrix) plane (Fig. 11), the hole mass remains con- 
stant with no angular dependence (See Table II). How- 
ever, as a consequence of a finite angular variation of the 
Fermi energy, which is pulled by the change in the den- 
sity of electron-like carriers, the hole lines present a small 
angular modulation. As seen in the figure, this theoret- 
ically expected angular variation is resolved by the ex- 
periment. Note also the doubling of the experimentally 
resolved hole peaks, which is a consequence of a slight 
misalignment. 

As far as the electron lines are concerned, the agree- 
ment between theory and experiment is less satisfactory. 
Note that the three distinct electron pockets become de- 
generate only when the field is along the trigonal axis. We 



found that the angular evolution of those Nernst peaks 
which can be attributed to electrons is in rough agree- 
ment with theory. However, this comparison is compli- 
cated by the persist ence of additional Nernst peaks of 
unknown origin^^'^fil^ which are hard to distinguish from 
the electron peaks. In the three Figures |9|10|11[ a spe- 
cific symbol marks those additional peaks. This symbol 
is used only when the resolved peak is clearly not one of 
the expected electron peaks. However, at this stage, in 
this field range and in particular when the field is almost 
oriented along the trigonal axis, we cannot definitely dis- 
tinguish between Nernst peaks associated with the elec- 
tron pockets and those which are unexpected. As the 
field increases and the number of expected electron lines 
decreases, it becomes more straightforward to clearly dis- 
tinguish between theoretically expected phase diagram 
and additional lines resolved by experiment. This was 
the case in a recent angle-dependent Nernst study ex- 
tended to high-magnetic fields by our groupP^. 

We have compared different sets of parameters in the 
model and opted for the set yielding the best agreement 
with experiment. These parameters yield an angular de- 
pendence for the electron lines, which still differ from 
the experimental results by a visible margin. This can 
be seen in Fig. |12[ which is a zoom on the region in 
the vicinity of the trigonal axis. The superposition be- 
tween the theoretical solid lines and the experimental 
color map is good but not perfect. Any attempt to re- 
pair this imperfection within the model used here led to a 
departure from the quite satisfactory agreement between 
theory and experiment on the subject of the angular vari- 
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FIG. 11: The Landau spectrum for a magnetic field rotating 
in the (binary, bisectrix) plane according to theory (fines) and 
experiment (symbols). Additional peaks refer to those which 
can be clearly identified as not associated with an electron 
Landau sub-level. 



ation of the hole lines. 

The imperfect agreement between theory and exper- 
iment regarding the electron lines may point to some 
deeper physics associated with electron interaction. A re- 
cent theory^ suggests that Coulomb interaction in pres- 
ence of strong mass anisotropy can lead to an inequality 
in the occupation of the three equivalent electron val- 
leys. Here, the relative agreement between theory and 
experiment puts an upper limit on any hypothetical in- 
equality in the occupancy of the three equivalent electron 
valleys. Given the impressive agreement between theory 
and experiment in the case of holes, however, the rela- 
tive disagreement in case of electrons (in particular for an 
arbitrarily-oriented magnetic field) may be speculatively 
attributed to a mild "spontaneous valley polarization" . 



VII. THE ZEEMAN SPLITTING OF HOLES 

A particularly important issue is the Zeeman splitting 
of the holes. When the field is along the trigonal axis, at 
low enough temperature, hole peaks are twin peaks, save 
for the one at the highest field at 9 T. This is a conse- 
quence of Zeeman splitting. The absence of splitting for 
the 9T peak associated with the n = 1 Landau level can 
be naturally understood if one assumes that the Landau 
indexes of the two adjacent peaks differ by two. In other 
words, the Zeeman energy creates a quasi-degeneracy be- 
tween the {n,+) and (n — 2,— ) Landau levels, when the 
field is parallel to the trigonal axis. On the other hand, 
when the field is perpendicular to the trigonal axis, the 




FIG. 12: Zoom on Color map of Sxy for a magnetic field near 
the trigonal axis. The magnetic field rotates either in the 
(trigonal, binary) plane {(f} — 0) (a) or in the plan (trigonal, 
bisectrix) plane {(j> = 90°). On top of the this color map, the 
theoretical Landau spectrum for the three electron pockets, ei 
(dotted line), e2 (solid line) and 63 (dash-dot line) are shown. 



Zeeman energy is vanishingly small. 

These features have been known since Smith and co- 
workergSl. From the position of the Zeeman-split peaks 
in the B || trigonal configuration, they concluded that 
Ez/fiiOc is cither 2.16 or 1.84. Indeed, the hole spec- 
trum for B//trigonal does not allow to distinguish be- 
tween these two possibilities. The Zeeman energy can be 
slightly larger or slightly lower than twice the cyclotron 
energy according to the way one indexes the twin peaks. 
Previous studies, focused on the spectrum in the vicin- 
ity of the trigonal axis, used either the former value^^ES 
or the latter^. On the other hand, there is no ambi- 
guity in indexing the Landau levels in the B _L trigonal 
configuration and, thus, the angular map allows to defi- 
nitely settle this issue. The model used here assumes that 
2.12 and, as seen in Fig. [9] and Fig. 10 



Ez/huJc _ 

cessfully reproduces the complex experimentally-resolved 
angular dependence of the hole lines and in particular the 
angle at which they cross each other. 

When the field becomes perpendicular to the trigonal 
axis the splitting is below the threshold of experimental 
resolution even in the subkelvin temperature range and 
in a field of 9 T. This can be seen in FigjTS] which tracks 
the evolution of the hole Zeeman splitting in the vicinity 
oi 9 ~ 90°. In other words, Ez/huJc is lower than the 
experimental resolution in agreement with theory (See 
table II). 



12 



= 87.7° 




92.7° 



FIG. 13: Nernst voltage as a function of B~^ for (f) = 90° as 
6 is scanned in the vicinity of 90°, which correspond to the 
bisectrix direction. The curve are shifted for clarity. The black 
arrows indicate the position of two hole peaks with opposite 
spins(4^ and 4~). Near 6 — 90°, experiment cannot resolve 
two distinct peaks. In other words, the Zeeman splitting be- 
comes smaller than the width of a peak. 



VIII. SUMMARY 

Our angle-dependent Nernst measurements map the 
Landau spectrum of bismuth up to 12 T. The experi- 
mental results are confronted with a theoretical model 
treating holes as quasi-particles with a parabolic disper- 
sion and electrons as Dirac fermions. Within this model, 
we found a set of parameters, which produce theoreti- 
cal results in good agreement with the experimentally 
resolved spectrum in the whole solid angle. The model: 

1. Fits almost the entire hole spectrum very well. 

2. Finds Ez/htUc — 2.12 for holes when the field is 
along the trigo nal ax is in agreement with the ex- 
perimental value^2I211_ 



3. Yields a satisfactory description of the level cross 
between 0^ and 1~ when the field rotates in the 
(Binary, Bisectrix), plane. 

4. Assumes a tilt angle for the electron pockets, 
6.2°, which is in agreement with previous low-field 
studie^i^. 

5. Concludes that Ij^" and 1^ meet at i? = 9T when 
the field is oriented along the trigonal axis in agree- 
ment with experiment. 

On the other hand, the model is not entirely success- 
ful. In particular, the fine structure near the trigonal 
axis according to theory and experiment are somewhat 
different. 

An important conclusion of this analysis is that when 
a field exceeding 2.5 T is applied along the bisectrix axis, 
all electrons reside in their lowest (j = 0) spin-polarized 
Landau level. 

Finally, let us add that additional Nernst peaks, un- 
expected in the non-interacting picture, were resolved on 
top of this complex spectrum. The explanation of their 
origin is a subject of ongoing research. 
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